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SUMMARY 


We use the correspondence between iterated integrals and a noncommuta- 
tive algebra to recast the given dynamical system from the time domain 
to the Laplace-Borel transform domain. Then, we show that the following 
algebraic criterion has to be satisfied for the onset of chaos: 


>25 <E P3)S u S = 0 

x„-»oo fc = o "'• x 0 


where Q is the generating power series of the trajectories, the symbol II 
is the shuffle product (le melange) of the noncommutative algebra, xq is 
a noncommutative variable, r is the correlation parameter. In the given 
equation, symbolic forms for both § and II can be obtained by use of 
one of the currently available symbolic languages such as PLI, REDUCE, 
MACSYMA. Hence, the criterion is a computer-algebraic one. 


INTRODUCTION 

Chaotic regimes of nonlinear dynamical systems are becoming more and 
more important as we find new applications. When we are interested in 
designing the operation of a dynamical system so as to avoid its chaotic 
regimes, we must know in advance the ranges of parameters for which the 
system’s response will be chaotic. This explains the need for a criterion. 

A dynamical system will respond according to the values of the physical 
parameters in play. We can always specify a space of physical parameters 
over which the system operates and imagine traversing paths within such 
a space. On any segment of a path, the system's dynamical response will 
be characterized by the long-term behavior of its trajectories in a phase 
space, and this characterization may change at critical points on the path. 
For example, over a certain segment of parameter space, the system’s re- 
sponse may become steady after a long time has elapsed. Corresponding 
trajectories in a phase space are attracted to a fixed point; accordingly, we 
call it a fixed-point attractor. On the other hand, over another segment of 
parameter space, the system’s long-term behavior may be periodic in time; 
corresponding trajectories in a phase space will approach a closed curve or 
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limit cycle. Since this curve is defined by one dimension, its length, we 
call it a one-dimensional attractor. Over yet another segment of parameter 
space, the system’s long-term response in time may contain two incommen- 
surate periods. Corresponding trajectories in phase space will be attracted 
to a two-dimensional figure, a torus; we call it a two-torus. Finally, one or 
more segments of parameter space may exist where the system’s long-term 
response appears to be random or chaotic. Corresponding trajectories in 
phase space will be attracted to a limiting set, the dimension of which is 
noninteger. We call this attracting set a strange attractor. To avoid chaotic 
regimes of the system, we must know the ranges of parameters over which 
the system’s response is characterizable by a strange attractor. 

We seek the necessary and sufficient conditions for the existence of a 
point in parameter space defining the onset of a chaotic regime in a given 
dynamical system. At the moment, we have two leads suggesting how we 
might arrive at such a criterion. First, if we knew how to tell when the 
sign of the largest Lyapunov coefficient became positive, we would have a 
criterion in our hands. At present, however, we do not know how to do 
this without computing the coefficients. The trouble here is that we must 
repeat the computation for each possible set of physical parameters, and we 
may never exhaust the possibilities. Numerical computations of this nature, 
therefore, are not suitable for developing a criterion which must be capable 
of picking out the relevant set of parameter values from the infinity of 
possible combinations. The second lead suggests itself from the observation 
that, in phase space, the trajectories will have zero autocorrelation at the 
onset of chaos. Here, rather than carry out numerical computations it is 
possible to derive an analytical expression for the autocorrelation function 
at the onset of chaos, and seek a criterion for the onset from the expression 
itself. To this end, it is convenient to transform from the time domain to 
the Laplace- Borel (LB) transform domain. 

Laplace-Borel transforms lend themselves to the treatment of iterated 
integrals. The connection between iterated integrals and nonlinear dynami- 
cal systems can be brought to light by giving the example of-let us say-the 
Duffing equation. We can integrate this equation twice and express the 
same amount of dynamics in the form of (twice-) repeated integrals. Re- 
peated integrals have been studied extensively by Chen in references 1, 2, 
and 3. There exists a formal correspondence between the repeated integrals 
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and a noncommutative algebra, as explained by Fliess and his colleagues in 
references 4 and 5. We shall use this correspondence to recast the dynam- 
ical system from its repeated-integral form to an algebraic form in the £8 
transform domain. The resulting set of algebraic equations can be solved, 
in principle, precisely for the physical parameters existant at the onset of 
a chaotic regime. 

This paper contains a derivation of the criterion along with a brief 
summary of the necessary mathematical background. Subsequent work 
should explore applications of the criterion to widely known dynamical 
systems as well as the novel flight dynamical system that we have discussed 
in references 6 through 9. 

This work was done under the NASA Ames Research Center/University 
of Santa Clara consortium, NCA2-48. 

BACKGROUND 

The following sections contain a brief summary on the use of iterated in- 
tegrals, a noncommutative algebra, and £8 transforms. More complete 
information can be found in references 1 through 5. 

Iterated Integrals 

Adopting the notation introduced by Chen (refs. 1, 2, and 3), we define 
an iterated integral as follows: 

l d ^"- d ^0 = f o J o d £jn—l • • • d £jo ( 1 ) 

and 

/V, = e,w - «,(») (2) 

JO 

where j = 0, 1, . . . , n and 0 < t < T and fo, £i, . . . , : [0, T] — > R are 

(n+ 1) continuous functions with bounded variations. We can associate the 
forcing terms of the dynamical system with iterated integrals as follows: Let 
it*; ui, U 2 , . . . u„ : [0, T] — > R; then a forcing term u, can be represented by 
the integral 

(3) 

J 0 
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where 


( 4 ) 

( 5 ) 


to{ T ) = T 

&(r) = / Ui{a)do 
Jo 

and i = 1, . . . , n. 


A Noncommutative Algebra 

We shall define a noncommutative algebra and its correspondence with the 
iterated integrals. Let X = {x 0 , . . . ,x n } be a finite set called an alphabet 
and let X = {u; 0 ,- • • , Wn} be a free monoid, another set, generated by X. 
If Wi is an element of X (w, G X), then w t is a word. The word w, consists 
of a finite sequence x, M . . . x, 0 of letters of an alphabet. 

The product of two words will be defined as a concatenation or juxtapo- 
sition of these two words. Concatenation is a noncommutative operation. 

The neutral element, which can be interpreted as an empty word, will 
be denoted by 1. 

The length of a word, say to*, will be denoted by | W( |, u>, G X, and it will 
be defined by the number of its letters. For example the word x 0 xixixix 0 
will have a length of 5; i.e., | x 0 xfx 0 |= 5 . The length of the neutral element 
is | 1 |= 0. 

The product in this algebra is called the shuffle product (le melange). 
We denote it by II and define it by induction on the length of words as: 


1 U 1 = 1 (6) 

xU1 = 11Ix = xVxGX (7) 

(xw) II ( x'w ') = x\w II (x'tt/)] + x'[(xu>) II w'] (8) 

Vx,x' G X and Vw,w' G X (9) 


We shall give an example of the shuffle product to familiarize the reader 
with this form. Let w\ = xo^i be the first word and W 2 = XiXq the second 
word. Let us determine tui II (read wi shuffled into w 2 ) : 

XqXj II XjXq = x 0 {xi II X 1 X 0 ) + Xi(x 0 Xi II x 0 ) (10) 
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The first term of the right-hand side generates the following terms: 

( 11 ) 

The second term of equation (11) generates the following terms: 

Xi(l II x 0 ) -I- x 0 (x! II 1) (12) 

Similarly, the second term of the right-hand side in equation (10) will give 
the following expression : 

x 0 (Xj II x 0 ) -I- X 0 (x 0 X! II 1) (13) 

which finally yields 


x 0 (xi(l II x 0 )) + x 0 (xi II 1) -(- x 0 x 0 xi (14) 

Next, we add like terms on the right hand side to obtain : 

X 0 (xiXiX 0 + X1X1X0 + XiXoXi) + Xi(x 0 XiX 0 + X 0 X 0 Xi -I- X 0 X 0 Xi) (15) 
Further regrouping yields : 


2x 0 XiXiX 0 + XQXiXoXi + XxXoXiXo + 2 X!XoX 0 Xi (16) 

Finally we obtain for the right-hand side of equation (10): 

2x 0 xjx 0 + x 0 xix 0 xi + xix 0 xix 0 + 2 xiXqXx (17) 

The reader should not be alarmed by the tediousness of this operation since 
a computer programmed for a symbolic language can carry out the shuffle 
product for us. 

Next we shall define the real algebra of formal polynomials and the real 
algebra of power series with real coefficients and noncommutative variables, 
Xj G X. Let 9?{x) be the 9?-algebraof a formal polynomial and 9? <C x > be 
the 9?-algebra of a power series with real coefficients and noncommutative 
variables. Let Q € 9? <C x » be an element of 9? <C x », which can be 
written as 

5 = £{[(5>w)]w I w € X) (18) 
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where (§,w) 6 3?. We define the following three operations ( addition, 
Cauchy multiplication and shuffle product ) on two elements of 3?, namely, 
: 


9i + 5i = + ( 5 2,w)]w 1 

w ex} 

(19) 

9i9i = Yl{\ 5Z (si,wi)(s2ie 2 )]iu 

| w e x} 

(20) 

ntj W2=W 

^1152 = X^{( S l’ U 'l)( 5 2,^2)«’l U W 2 1 

wi,w 2 e x} 

(21) 


Let us give an interpretation of the meaning of a noncommutative power 
series § which is an element of 3? <C i >. This noncommutative power 
series actually defines a causal (nonanticipative) functional of the forcing 
functions u, of the dynamical system if we replace the word . . . ij 0 by 
the corresponding iterated integral / 0 ‘ . . . d£y 0 - Let y(t : u l5 . . . , u n ) be a 

causal functional. We can write this causal functional as 

" f* 

J/(f i ^i> * • • i ^n) (9 ■> 1 ) T ^ ^ 'y ' {9 •> ■ • • Zjo) I ■ • • d£jo (22) 

H>0 j0,...,jn=0 ® 

Such a causal functional is said to be analytic with generating power series 
9- Equation (22) should be considered as a representation of a causal func- 
tional in terms of a noncommutative generating power series and iterated 
integrals. 

Let us consider an example to better understand what all this means. 
Let u\ represent the forcing function of a dynamical system. Let u 3 be 
given as 


Uj = t° 

(23) 

Then we shall have (eqs. 3-5) 


[ d£ 0 d£ 1 = f dr [ o a do 
Jo Jo Jo 

(24) 

= r r<lTl dr 
J 0 a + 1 

(25) 

t a + 2 

(a + l)(a + 2) 

(26) 
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Now let us consider the following integral instead 


/ dfidfo = f o a do [ dr = [ o a .odo 
'o Jo Jo Jo 

(27) 

= f o a+x do 

Jo 

(28) 

t<*+2 

a + 2 

(29) 


These two repeated integrals are different if and only if (iff) a ^ 0 in the 
same way that the corresponding series xoXj ^ XiXo- This is what we mean 
by noncommutativity. 

We shall give some fundamental theorems and remarks from references 
4 and 5 and urge the reader to consult the references for the details of the 
proofs: 

1. A causal functional is said to be analytic if it is defined by a noncom- 
mutative formal power series called the generating power series. 

2. The notion of an analytic causal functional generalizes in some sense 
the notion of an analytic function. We can represent an analytic 
function by its Taylor-series expansion and in a similar fashion we 
can represent an analytic functional by its generating power series. 

3. The product of two analytic causal functionals is a functional of the 
same kind, the generating power series of which is the shuffle product 
of the two generating power series. 

Laplace-Borel Transforms 

In this section we shall define the CB transformation. Let us consider an 


analytic function h(t): 

h(t) = '£h n t - 

n> 0 7t * 

(30) 

We also have 

t n 

n! 



/ dr n I dr n -i . .. f dr , 
Jo J 0 Jo 

(31) 
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Hence, we can rewrite h(t) as 

h{t) = Yl k rtfdT n f ... f dr i (32) 

r^O J ° J ° J ° 

Let us define the generating power series of noncommutative variable x 0 as 

$ = !>»*; o ( 33 ) 

n>0 

We say that the ZB transformation of h(t) is Q and we write 

£B(h{ t)) = S='£h„x" 0 (34) 

n>0 

We have introduced a number of new concepts (new even to some math- 
ematicians). An example may help to fix and identify them. Let us consider 
the Duffing equation: 

y(t) + ay{t) + y(t) + 0y 3 (t) = Uj(t) (35) 

The nonlinearity is a cubic polynomial, so that the equation is sometimes 
referred to as a cubic harmonic oscillator. To obtain the repeated integral 
form, we integrate twice : 

1 . Integrate from 0 to t 

ylo + <*ylo +/ y(r)rfr -H A y{r) 3 dr = f u^dr (36) 
Jo Jo Jo 

and rearrange the terms to 

y(t) + ay(t ) + f y{r)dT + /? f y(r) 3 dr = y(0) + ay( 0) + f «i(r)dr 
Jo Jo Jo 

(37) 


2. Integrate again 

y[t)-ha [ y(r)dT + [ dr 2 f y(r 1 )dr 1 + 0 f dr 2 f y^) 3 ^, = (38) 

Jo Jo Jo Jo Jo 

/ y(0)dr + f a(y[0))dT + y(0) + f dr 2 f Ui(T 1 )dr 1 (39) 
Jo Jo Jo Jo 
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3. Use the correspondence between the iterated integrals and noncom- 
mutative variables to obtain 

9 + ax 0 9 + *19 + 0*1 5 II5U^ = &+(o:a4- b)x 0 + x 0 x 1 (40) 
where 

a = y(0) 

6 = y(0) 

are the initial conditions. We can rewrite this as 

(1 + ax o + x 2 0 )9 - b + (aa + b)x 0 + x 0 Xi - 0xl9 II 9 U 9 (43) 

The important point is that in the ZB transform domain, the given non- 
linear dynamical system is now represented as an algebraic equation. This 
algebraic equation can be solved, for example, by iteration. We consider it 
very important to acknowledge that, indeed, we are able to reduce a given 
nonlinear differential equation to an algebraic equation in much the same 
way that we use ordinary integral transforms (Laplace or Fourier) to reduce 
a given linear differential equation to an algebraic equation. An iterative 
scheme to solve equation (43) can be formulated as follows: Let us try a 
solution of the form 


(41) 

(42) 


9 — 9o + 09i + ft 2 9% + • • • 

Substitution into the original equation yields 

6 -f- (aa + 6)xq 4" XqXj 


9o - 


1 + ax 0 + Xq 


$i = - 


1 + ax 0 + xg 


(£oU£olI£o) 


(44) 

(45) 

(46) 


r - x o 

t (2 

1 + axo + Xq 
We can rewrite 9i 35 


(9o u 9i n §o + 9i u 9o u 9o + 9o u 9o u 9i) (47) 


9 2 - - 


1 -f ax 0 + Xq 


(3^oU5oII^i) 


(48) 
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1 + ax 0 + Xq 

x {9oU9 2 ugo+$iugoU$o+92ug 0 u$o+$oug i ug 1 +g 1 u9 0 ug 1 ) 

( 49 ) 


Let us introduce 


Then we have 


S (x 0 ) = 


1 + ai 0 + x 


2 

o 


£o = -S(xo)(- + 

Xq 


aa + b 

Xq 



5\ — S (xojpo U So U 9o 
9i ~ 25(xo)^o u So if 5\ 


( 50 ) 

( 51 ) 

( 52 ) 

( 53 ) 


5s = S(*o)(5oU5s1J5o~5iU5oU5o+5 2 U5oU5c,+ SoU5iU5 1 +.5 1 U5 0 115 1 ) 

( 54 ) 

We cannot ignore the connection between these iterations and Feynman 
diagrams, which enable one to figure out the next iteration if one does not 
want to use a computer. We repeat that a computer can be programmed 
to carry out such iterations for us up to whatever order we request. 

To return to the time domain we must: (a) decompose Si into partial 
fractions and (b) use the following lemma obtained in reference 1. 

Lemma: The rational fraction (l — ax 0 )~ p corresponds to the exponential 
polynomial 

J \ aW 

p - 1 / y ! 

A List of LB Transforms 

We give a partial list of LB transforms: 



( 55 ) 


CB(unitstep) = 1 


CB{e at ) = — - — 
1 — ax o 
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1 


£B{e at {l+at)) = 


(1 - ax o) 2 


£B|e"(l + al+?^)] = (l-az„) s 


mz (— ) ^)e of ] = (1 - a 0 z 0 r 

t=0 v 1 / *• 


£B($tnwt) = - 

1 + W 2 X o 

£S(C0 ~ ,) = FT ^ 


£B(sinhwt) 


rvxp 

1 — W 2 Xq 


£B(coshwt) = 


CHAOTIC REGIMES OF NONLINEAR 
DYNAMICAL SYSTEMS 

Chaotic regimes correspond to strange attractors. Strange attractors are 
limiting sets with noninteger dimension in the phase space of the trajec- 
tories of the dynamical system. Numerically speaking, given the physical 
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parameters we can compute: (a) the power spectrum; (b) the Lyapunov co- 
efficients; or (c) the dimension of the attractor to check whether the regime 
is chaotic or not. For the random-like response which is characteristic of the 
so-called chaotic regime, the computation (a) would result in a broad band, 
(b) would result in a positive Lyapunov coefficient, and (c) would result in 
a noninteger dimension. A fourth check would be (d) the autocorrelation 
function AC(t). The autocorrelation function AC(t) vanishes in the limit 
as t - — » oc for chaotic regimes. Trajectories in the phase space that start 
out infinitesimally close ultimately diverge exponentially, and hence fail to 
correlate at large times. More precisely, the autocorrelation is a measure 
of the similarity of a trajectory at a given time t with its value at a later 
time t + t. We define AC(t) as the arithmetic mean of a large number of 
products such as x(t).x(t + r); i.e., 

AC(t) = — - — f x(t).x(t + r)dt (56) 

t<i — 1 1 Jh 

In a more compact way we can write 

AC(t) = (x(t).x(t + r)) (57) 

We call AC(t) a temporal autocorrelation function. By varying the inter- 
val r, we can construct AC(t). It defines the degree of similarity of the 
trajectory x(t) with itself as time evolves. 

The Wiener-Kintchine theorem states that AC(r) is the Fourier trans- 
form of the power spectrum. As a result of this we see that for the regimes 
represented by an attractor which is: (a) a fixed point; (b) a limit cycle; 
and (c) an n-torus we have 


lim AC{t) ? 0 (58) 

T—*00 

since in these three cases the power spectrum is formed of distinct rays. 
In other words, periodic or quasiperiodic trajectories keep their internal 
similarity with the evolution of time. This means that the behavior of the 
system is predictable. In contrast, for a chaotic regime, where the power 
spectrum has a broad band (fig. 1) 

lim AC[t ) = 0 (59) 

T — ► OO v ' 
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As figure 1 illustrates, the temporal similarity of the trajectory with 
itself tends towards zero at sufficiently large times. It follows that know- 
ing x(t) for a long lapse of time (as long as one wishes) does not allow 
the prediction of the future of x(<) . We can say that the chaotic regime is 
unpredictable because of the progressive loss of internal similarity of tra- 
jectories. We observe a loss of memory in terms of initial conditions (7C) 
with attractors of the fixed point, limit cycle, and n-tori types since for 
each of these types, trajectories originating from diverse initial conditions 
outside the attractor converge on the attractor. For example, any initial 
condition will result in the same constant steady motion if the attractor is 
a fixed point. In chaotic regimes, we have just the opposite: almost the 
same initial conditions result in different final states. Figure 2 depicts the 
phase space in which two trajectories that are initially close diverge from 
each other and hence lose their similarity in a finite time. 

We shall give the following definition of chaos and its connection with 
the property called sensitivity to initial conditions (SIC) : 

1. Definition: If a regime is represented by an attractor such that in- 
finitesimally close neighboring trajectories diverge exponentially, such 
a regime is characterized as chaotic. The sensitivity to initial condi- 
tions implies an exponential amplification of errors or uncertainties 
in the IC and it characterizes the regime as chaotic. 

To summarize, the existence of a chaotic regime of a dissipative nonlinear 
dynamical system is characterized by any one of the following conditions: 
(a) an attractor with a noninteger dimension; (b) a sensitivity to initial 
conditions; (c) a broad-band spectrum; and (d) a vanishing autocorrelation 
function. 

We shall make use of the last condition for our criterion. We have 
observed that the ZB transformation will enable us to have a semianalytic 
criterion. The following basic properties of ZB transforms are needed in 
the formulation of the criterion: 

Shifting Theorem: The ZB transformation of a function whose argu- 
ment is shifted by an amount r is related to its ZB transform as follows; 

«!/<* + ')] -(Ej^W<01 («) 
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Proof: Let us write the Taylor-series expansion for this function; 

/(« + r) = /(0 + <7*M + ^/"(<) + ~- (61) 

Since we started with a function for which there exists an LB transforma- 
tion, we can write 

m = E <■-£ (62) 

n> 0 n ' 

or, there is a corresponding generating power series § of the noncommuta- 
tive variable x 0 ; 

P = (63) 

n>0 

This C) by definition is the LB transformation of the original function, f(t). 
The LB transformation of the derivative of f(t) follows immediately as; 


/'(<) = E 

n>0 n ' 

(64) 


(65) 

/(0_ §“”(»- i)’ 

(66) 

Hence its generating power series will be 


E 

n>0 

(67) 

or 

X ° n>0 

(68) 

or we write 

Xq 

(69) 

Now we can write down the LB transformation of 


f(t + t) = f(t) + t f 1 (t) -|- —/"(«) 

(70) 
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£B[f{t + r)] = Y1 Q n x 0 + — E a « x o + ^iE a " x 0 + • • • ( 71 ) 

n >0 X ° n > 0 Z ' X 0 n > 0 

£B[/(* + r)] = fl + — + — 2 + • • •) EZ a " x o ( 72 ) 

V x o * lx o / „>o 

£B[/(t + r)]=(5:-^)5:a n x" (73) 

V*>0 k - X 0 J n> 0 

or we write 

£«|/(i+r)]-fe^)"l/WI (74) 

which completes the proof of the shifting theorem. 

Next we show a limiting property under the name of the limiting theo- 
rem. 

Limiting Theorem: The temporal limit as t — > oo corresponds to the 
limit of the noncommutative variable of the L B transformation as x 0 — * 
oo. 

Proof: The proof follows from the definition of the LB transform. If we 
have an analytic function 

/(«) = £ «»£ ( 75 ) 

n >0 U - 

then its LB transform is given by 

£»[/(<)]=£<*»*?• ( 76 ) 

n>0 

Therefore, if t — ► oo then in the transform domain x 0 — *• oo. With the 
help of these theorems we are ready to present the main theorem of this 
paper, which represents the criterion for the onset of chaos. 

Main Theorem: The following algebraic criterion has to be satisfied for 
the onset of chaos: 

Jto (i:^Wu£ = ° (77) 

z„->oo \k > 0 K ' x 0 J 

where Q is the generating power series for the trajectories of the nonlinear 
dynamical system. In other words, Q is the LB transformation of the 
trajectories of the given nonlinear dynamical system. 
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Proof: We adopt as definition of the autocorrelation function 


AC(t ) = ^ J ^x{t).x{t + r)dt (78) 

As we discussed previously, we must have 

lim AC(t ) = 0 

T — >OC 


for the onset of chaos. 

Let us take the £B transformation of both sides, 

£B\AC(t ) ] = ^lim^xo-CS [x(0] U £B [z(f + r)] (79) 

Notice that we have used a result in reference 1 to take the £B transforma- 
tion of a product of two functions. The £B transformation of the product 
is equal to the shuffle product of the individual transformations. Let us 
denote £B [x(f)J by Q. Then by the shifting theorem we have 

" w ‘ + r,l = fe6)* 

Hence the criterion for the onset of chaos becomes 

lim |£iri)s u S = 0 
which completes the proof of the main theorem. 

CONCLUDING REMARKS 

The criterion obtained as a main theorem is a computer-algebraic one. 
It is noteworthy that the algebra required in obtaining the two generating 
power series and their shuffle product, together with the two limits, all can 
be done on a computer with a symbolic language. At the moment we have 
the option of using any of the following symbolic languages, PLI, REDUCE, 
MACSYMA, or LISP. 


(80) 


(81) 
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Although the criterion is simple in concept, its application by no means 
will be a trivial affair. A considerable effort doubtless will be necessary to 
handle the difficulty of bifurcation points. These are points in parameter 
space where analyticity will be lost, and hence where the generalized series 
expansions will not be valid. Such points will reveal themselves beforehand 
by slowing down the rate of convergence of the iterations. We anticipate 
the necessity of considering slow convergence as a signal to change the set 
of parameters so as to skip over a bifurcation point. Since we are not 
interested in the identification of such a point, it should not be necessary 
to come too close to it. We assume that there will be a finite and a small 
number of bifurcation points which it should be possible to skip successively 
if necessary, bringing us ultimately to a parameter regime that borders on 
the object of our search. This is one of the difficulties we can foresee. There 
may well be others. 

Our criterion, which is computer-algebraic, can be used to characterize 
the ranges of the physical parameters for which chaotic regimes will take 
place. Hence, it can be used in designing the operation of a nonlinear 
dynamical system to avoid such regimes. In other words, we can control 
chaos in our nonlinear dynamical system by staying outside of the ranges 
defined by our criterion. 
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Figure 2.- Two initially close trajectories. 
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